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(57) Flowmeter for the measurement of a flow char- 
acteristic of a fluid using a signal generator and at least 
one signal sensor for sensing a signal as it is carried in 


FIG. I 



the flow direction of the fluid in which signal generation 
and processing are operated to provide accurate and 
rapid flow measurements over a wide dynamic flow 
range for diverse fluid compositions. 
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Description 

The present invention is generally directed to a flow- 
meter for the measurement of a flow characteristic of a 
fluid using a signal generator and at least one signal 
sensor for generating and sensing a signal as it is car- 
ried by the fluid in the direction of flow 

Tlme-of-Flight (TOF) measurement is a flow sens- 
ing technique in which the flow of a fluid such as a gas 
or liquid is measured by timing the passage of a gener- 
ated signal, such as a thermal pulse, as it is carried be- 
tween two points over a known distance in the direction 
of the fluid flow The signal is convected by the flow of 
the fluid and the time it takes to pass between the two 
points is a measure of a flow characteristic of the fluid. 
The time of transit or flight of the signal is an indication 
of the mean velocity between the two measurement 
points which allows the volumetric flowrate of the fluid 
as well as other flow characteristics to be determined. 

The signal that may be used can be generated from 
a number of sources including heaters, acoustic gener- 
ators, optical generators and the like. What is required 
is that the signal be conveyed downstream by the fluid 
flow so that it can be detected by a sensing device. 

Thermal time-of -flight flowmeters have been known 
since 1948. LS.G. Kovasznay, "Hotwire Investigation of 
the Wake Behind Cylinders at Low Reynolds Numbers", 
Proceeding Royal Society vol. 198:174-190(1948) pro- 
vides a time-of-flight flow measurement for the flow of 
gas in a wind tunnel. A stationary heater is positioned 
normal to the mean flow upstream of a movable sensor. 
The spatial wavelength of the convected heat pattern is 
measured by moving a hot wire anemometer in the flow 
while knowing the heater driving frequency and compar- 
ing the sensor signal to the same. From this, the mean 
convection velocity may be calculated. 

Andrew B. Bauer, "Direct Measurement of Velocity 
by Hot-wire Anemometry, AIAA Journal, 3(6): 1 1 89-1 1 91 
(June 1 965) describes a heater-sensor probe arrange- 
ment for taking gas velocity measurements. The device 
is intended for very high velocity flow measurements un- 
der turbulent conditions. The heater and sensor com- 
prise short lengths of very fine wires made of platinum 
and rhodium which are parallel to one another and nor- 
mal to the direction of flow of the fluid. 

LJ.S. Bradbury et al., "A pulsed-wire Technique for 
Velocity Measurements in Highly Turbulent Flows", 
Journal of Fluid Mechanics vol. 49 (part 4) pp. 657-691 
(1 971 ) discloses a probe in which the sensor was posi- 
tioned at right angles to the heater. The authors discuss 
the effects of thermal diffusion and errors that may arise 
as a result, particularly at low flow rates and sensor sig- 
nal pattern stability which is responsible for errors at 
high flow rates. 

Unlike a velocity probe, a volumetric flow device 
measures the flow of a fluid passing through a conduit 
as described in C. A. Mosse "The Design and Evaluation 
of Devices for Measuring Breathing Tidal and Minute 


Volumes, Masters Thesis Cranfietd Institute of Tech. 
(Sept 1980). Mosse employs principals utilised with ve- 
locity probes for the problem of volumetric flowrate. The 
volumetric flow sensor described by Mosse was poten- 

5 tiaily suitable for neonatal volumetric flowrate monitor- 
ing. To counteract some of the noise that interfered with 
the heater-sensor arrangement, Mosse added an extra 
sensor wire at the sensor axial location but displaced 
from the first sensor wire. The dual sensor system was 

10 merely used to collect noise signals from the flow and 
did not intercept the thermal pulse for the purpose of 
generating a timing signal. 

C.A. Mosse et al., "Microprocessor-based Time-of- 
Flight Respirometer", Med. Biol. Eng. Comp., Vol. 25: 

is 34-40 (1987) incorporated a microprocessor in the 
processing electronics for real time operation of a trans- 
ducer using a single heater and two sensor wires. The 
heater-sensor pair were parallel to each other and on 
the transducer axis and hence parallel to the direction 

20 of the mean fluid flow. 

Goldberg, U.S. Patent No. 4,938,079 discloses a 
thermal transit time flow measurement system for liq- 
uids and particularly to a device for IV infusion. The de- 
vice incorporates a microprocessor for controlling the 

25 heater and signal processing. The system employs at 
least one heater and at least one sensor with two and 
three sensors being specifically disclosed. The heater 
and the sensors are located external to the conduit. Mi- 
crowave energy is used to generate the thermal pulse. 

30 The power from the heater is applied for a specified pe- 
riod such that the liquid flow experiences a typical tem- 
perature rise of a few degrees Centigrade. The second 
sensor is thermally isolated from the heater and senses 
the resulting thermal pulse as it is carried by the liquid. 

35 The flow measurement involves comparing the signals 
from the first and second sensors as a function of time. 
The transit time measurement is then related to the flow- 
rate through the use of the microprocessor. 

S. A. Abesekera et al. "Liquid Flow Measurement by 

40 Cross Correlation of Temperature Fluctuations", Trans- 
actions of the Institute of Measurement and Control 5 
(11) pp.435-439 (Nov. 1972) discloses a cross correla- 
tion liquid flowmeter for the measurement of steady and 
pulsating, laminar and turbulent flows. 

45 Flowmeters of the type discussed above suffer f rom 
a number of disadvantages. In many cases, the flowm- 
eters are dependent on the type of fluid that passes 
through the conduit. For example, some flowmeters are 
only applicable to liquids or gases but not both. In addi- 

50 tion, such flowmeters are limited by the composition of 
the gas or liquid flowing therein. As a result, the flowm- 
eter must be recalibrated or another flowmeter em- 
ployed when there is a change in the fluid composition. 
Another disadvantage of prior art systems concerns the 

55 accuracy of the flow measurement. Such systems are 
inaccurate or are accurate only within a narrow range of 
flow rates. Thus, some prior art systems may be em- 
ployed only when the difference between the minimum 
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and maximum flows through the conduit is within a spec- 
ified narrow range. 

Another disadvantage of prior art systems is that 
they are generally adapted to measure flowrates only 
for turbulent flows or laminar flows but not both. A still 
further disadvantage of prior art systems is that they lack 
the speed of response necessary for applications such 
as healthcare devices that require very rapid flowrate 
measurements such as employed for the administration 
of anaesthetics. Furthermore, it is typical for prior art de- 
vices to measure flow in only a single direction (i.e. uni- 
directional). However, some applications such as in 
healthcare devices, require bi-directional flow measure- 
ments (i.e. the measurement of flowrates of fluids going 
into and out of a patient). 

A flowmeter useful in a commercial setting such as 
in the healthcare industry must provide information 
about the flow characteristics of a fluid that is reliable 
and timely. The performance of a flowmeter is depend- 
ent on a number of performance characteristics. Includ- 
ed among these are fluid composition independence 
(especially for gas flow measurements), speed of re- 
sponse, flow dynamic range, accuracy including repeat- 
ability and linearity, and minimisation of crosstalk inter- 
ference. Other performance characteristics which im- 
pact on the operation of a flowmeter include ruggedness 
of the flowmeter, compactness, and the power required 
to operate the device. 

As used herein the term fluid composition inde- 
pendence* shall generally mean that the device can pro- 
vide reliable and timely information of flow characteris- 
tics regardless of the composition of the fluid. In partic- 
ular, fluid composition independence provides a flowm- 
eter which performs in a desirable manner regardless 
in the change of density, thermal conductivity, viscosity 
and/or specific heat of the fluid which is to be measured. 

The term 'speed of response" shall generally mean 
the time it takes for the flowmeter to detect a change in 
flow rate. 

The term "accuracy" as used herein includes two 
components; namely, repeatability and linearity. The 
term "repeatability" shall generally mean herein that the 
flowmeter is capable of providing the same measure- 
ment for the same flow over a set period of time. In other 
words, if the flowrate remains constant, the readout of 
the flowmeter will remain constant and not fluctuate sig- 
nificantly (e.g. by more than 5%). 

The term "linearity" is a comparison between the 
mean flowrate at a given point in time and the actual 
flowrate. A plot of the mean flowrate versus the actual 
flowrate for a high performance flowmeter as described 
herein should result in essentially a straight line re- 
sponse. 

The term "flow dynamic range" generally means the 
range of flowrates that can be measured by the flowm- 
eter. A "wide flow dynamic range" as used herein in as- 
sociation with the present invention is one in which both 
very low flowrates (approaching zero) and very high flow 


rates can be measured by the same flowmeter. More 
specifically, a wide dynamic flow range generally in- 
cludes a flow range in which the ratio of the lowest to 
the highest flowrates is at least about 10:1, preferably 

s at least about 50:1. 

It would be a significant advance in the art of flow 
measurement to provide a flowmeter which can provide 
accurate and rapid flow measurements over a wide flow 
dynamic range for diverse fluid compositions in a cost 

10 effective and efficient manner. 

According to a first aspect of the present invention, 
a device for measuring a flow characteristic of a gas 
comprises: 

is (a) at least one signal generating means positioned 
in proximity to the flowpath of the gas for generating 
a detectable signal and for delivering the detectable 
signal within the flowpath of the gas; and 

20 (b) at least one pair of longitudinally spaced apart 
sensors positioned downstream of the signal gen- 
erating means for sensing the detectable signal as 
a measure of the time of flight of the detectable sig- 
nal between the pair of sensors. 

25 

The configuration of heater-sensor-sensor as em- 
ployed in the present invention has intrinsic advantages 
related to independence of the flow measurement from 
gas composition variations. 

30 in a preferred embodiment the signal generating 
means is adapted to generate a thermal pulse having a 
peak temperature in which the shape of the pulse is sub- 
stantially uniform about the peak temperature. By con- 
trolling the shape of the thermal pulse as described 

35 hereinafter, errors due to variations in gas compositions 
are reduced or eliminated. 

In a related aspect of the invention, there is provided 
at least one line sensor positioned off-axis of the flow- 
path which senses the signals along a line traversing 

40 the direction of the flow of the fluid. This aspect of the 
invention improves linearity and therefore accuracy and 
makes it possible for the present flowmeter to measure 
flow characteristics in both turbulent and laminar flow 
regimes. 

4S in one embodiment of the present invention, there 
is provided a flowmeter having a heater which can gen- 
erate thermal pulses at a rate where it is not necessary 
for the heater to cool down to the original starting tem- 
perature. This system enables the more rapid genera- 

so tion of thermal pulses and thereby faster response 
times. In a preferred aspect of the invention, the heater 
is a thin film conductor comprising a solid electrically 
non-conducting core and a thin metallic film covering at 
least a portion of the core. Similarly constructed sensors 

55 can also be employed. 

In a still further embodiment, the signal which is 
generated from the signal generating means comprises 
the simultaneous combination of at least one higher f re- 
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quency component and at least one tower frequency 
component. The employment of multiple frequencies 
provides a flowmeter capable of measuring fluids in a 
wide flow dynamic range (i.e. very low flowrates to very 
high flowrates) without time consuming calibrations or 5 
the use of multiple flowmeters. 

There is also provided in accordance with the 
present invention a system for processing simultaneous 
multiple frequency signals to determine the time delay 
difference between two instances of the signal both ac- 10 
curately and unambiguously over a wide range related 
to the periods of the frequency components. This sys- 
tem eliminates or at least reduces phase ambiguity er- 
rors particularly at low flowrates. 

In yet a further embodiment of the present invention is 
the heater generates alternating power pulses, prefera- 
bly such pulses having alternate polarities. The gener- 
ation of detectable signals in this manner shifts the elec- 
trical power spectral power to half the thermal pulse sig- 
nal frequency and thereby serves to eliminate crosstalk 20 
interference when the signal is processed. 

Related to this aspect of the invention, the crosstalk 
interference signal can be used to eliminate phase delay 
errors due to normal differences due to manufacturing 
tolerances in the signal processing components of each 25 
channel of the signal processing system. In particular, 
the signal generation means generates not only a signal 
to be detected but also inherently generates a crosstalk 
interference signal. Means of sensing the crosstalk in- 
terference signal and for employing the crosstalk inter- 30 
ference signal to eliminate any phase delay errors of this 
type in the detectable signal is also provided. 

The flowmeter of the present invention can also em- 
ploy systems for controlling the temperature of the heat- 
ers so as to (a) ensure that they do not exceed a prede- 3S 
termined maximum temperature and/or (b) are operated 
within a selected range of temperatures. Such systems 
are desirable for flowmeters which operate under a wide 
flow dynamic range such as in the healthcare industry 
in devices for measuring the inflow and outflow of gases 40 
by a patient. 

The accuracy of the present flowmeter can be en- 
hanced by having a multiplanar arrangement for the sig- 
nal generating and signal sensing devices. In particular, 
the flowmeter of the present invention can include a plu- *s 
ralrty of signal generating and sensing units wherein 
each unit comprises at least one signal generator and 
at least one signal sensor. Each of the units is positioned 
in a plane along the mean flow direction of the fluid. The 
planes are parallel to each other so that the units are so 
spaced apart and parallel to each other in the mean flow 
direction of the fluid. In addition, the signal generator 
and sensing devices of one unit are in the same longi- 
tudinal position of another unit. 

In a still further aspect of the present invention, the ss 
flowmeter can detect a flow characteristic of the fluid re- 
gardless of the direction that the fluid is flowing. For ex- 
ample, if the flowmeter is positioned within a conduit 


such as associated with an intubation tube, the flow di- 
rection of the fluid can be going into the patient or in the 
opposite direction, going out of the patient. In one aspect 
of the present invention, the flowrate of the fluid can be 
measured in either direction. This is accomplished by 
providing first and second signal generating means for 
generating a detectable signal within the fluid flow so 
that the same can be detected downstream by the sen- 
sors. 

In another aspect of the invention, there is provided 
a method of determining the time delay difference be- 
tween two instances of the same signal unambiguously, 
said signal containing at least two known frequency 
components. The method comprises cross correlating 
the signal over a correlation interval substantially equal 
to the period of the lowest frequency component and 
locating the peak of the cross correlation function to 
thereby obtain the time delay. This aspect of the inven- 
tion is an efficient cross correlation technique for deter- 
mining the thermal signal transit time unambiguously 
where the signal comprises a range of known simulta- 
neous frequencies in order to achieve wide dynamic 
range performance. The efficiency of the technique is 
such that it is helpful for smooth operation in an auto- 
matic real time device. 

The multiple aspects of the present invention en- 
hance one or more of the components of performance 
to obtain a flowmeter device which can be independent 
of fluid composition, provide for rapid speed of response 
over a wide flow dynamic range. In addition, the accu- 
racy of the measurement including repeatability and lin- 
earity are improved over systems known in the art. 

Embodiments of the invention will now be de- 
scribed, by way of example, reference being made to 
the Figures of the accompanying diagrammatic draw- 
ings in which: - 

Figure 1 is a schematic view of an embodiment of 
the flowmeter system of the present invention; 

Figure 2 is a cross-sectional view of an embodiment 
of a transducer which can be used in the flowmeter 
system of the present invention showing the relative 
position of the heaters and sensors; 

Figure 3 is a front view of the transducer shown in 
Figure 2; 

Figure 4 is a rear view of the transducer shown in 
Figure 2; 

Figure 5 is a cross-sectional view of a thin film con- 
ductor which may be employed as a heater or sen- 
sor in the present invention; 

Figure 6 is a schematic view of one embodiment of 
the invention for generating a thermal pulse com- 
prising a frequency modulated and/or amplitude 
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modulated waveform component; 

Figure 7 is a graph showing a signal spectrum gen- 
erated through the use of a frequency modulated 
heater voltage waveform; 

Figure 8 is a schematic view of one embodiment of 
the signal processing system of the present inven- 
tion; 

Figure 9A is a flow chart showing the flow of digital 
data involved in the processing of a detectable sig- 
nal from the sensor to a readout of information re- 
garding the flow characteristics of a fluid; 

Figure 9B is a graph showing a cross correlation 
function of a signal containing two different frequen- 
cies and a delayed version of that signal where the 
delay is greater than the period of the high frequen- 
cy component; 

Figure 10 is a flow chart showing the processing of 
the high and low signal strength measurements as 
shown in Figure 9A to determine flow direction, in- 
strument flow zero and the like; 

Figure 11 is a schematic view of a signal processing 
system using zero crossing timers; 

Figure 1 2 is a flowchart showing a master algorithm 
for processing the signal in the embodiment of Fig- 
ure 9A; 

Figure 13 is a schematic view of a still further em- 
bodiment of the invention including two systems for 
controlling the temperature of the heater; 

Figure 14 is a schematic view of a divider circuit 
used in the embodiment of Figure 1 3; 

Figure 1 5 is a graph showing the inputs and outputs 
associated with the timed divider voltage sampling 
of Figure 13; 

Figure 16 is a flow chart of an algorithm for opera- 
tion of the temperature heater control systems 
shown in Figure 1 3 as well as optional amplitude 
modulation generation; 

Figure 17 is a schematic view of an embodiment of 
the invention using a biplanar arrangement of heat- 
ers and sensors; and 

Figure 18 is a schematic view of a further embodi- 
ment of the signal processing system of the present 
invention particularly adapted for a biplanar ar- 
rangement of heaters and sensors. 


The present invention is directed to a flowmeter sys- 
tem in which the flow of a fluid, whether gas or liquid, 
can be measured to obtain desired flow characteristics 
including volumetric flowrate, or velocity. As a conse- 

5 quence of the present invention, measurements of 
these flow characteristics are obtained with improved 
performance over flowmeters known in the art. In par- 
ticular, the present invention provides a flowmeter sys- 
tem which can measure fluids regardless of composition 

10 over a wide flow dynamic range with exceptional speed 
of response and accuracy as defined herein. 

In general, the flowmeter system provides for the 
generation of a signal, which can be created, for exam- 
ple, by varying the temperature of a heater proximate to 

is the fluid and then relating the temperature variation of 
the fluid into a measurement of a flow characteristic of 
a fluid, such as can be performed with a heater-sensor 
system. The flowmeter system preferably includes a 
combination of a signal generator and a signal sensor, 

20 preferably a pair of downstream sensors for sensing the 
generated signal downstream of the fluid flow. A fre- 
quency encoding system is provided to generate detect- 
able signals under controlled conditions which can be 
readily sensed by the signal sensors. There is also pro- 

25 vided a signal processing system by which the signals 
sensed by the sensors are processed to obtain a flow 
characteristic of the fluid. 

The present flowmeter can sense any type of fluid 
including gases and liquids. While liquids can be sensed 

30 in accordance with the present invention, a description 
of the present system will be made with reference to 
gases only. It should be understood however that the 
system, unless otherwise specified, can be applied to 
liquids as well. 

35 The signal generating means employed in the 
present invention can include devices for generating a 
wide range of signals including thermal pulses, acous- 
tical signals, and the like. All that is required is that the 
signal generating means generate a signal which can 

40 be carried with the fluid downstream to a sensor. The 
signal sensing means can be any device which can 
sense or detect the signal generated by the signal sens- 
ing means which has been carried in the flow direction 
of the fluid. 

45 u will be understood that a variety of different signal 
sources can be used. By way of example and for con- 
venience only, the detailed description of the invention 
will refer to a preferred form of the invention including a 
signal sensing means (a heater) which is capable of 

50 generating a detectable signal (a thermal pulse) which 
is sensed by a sensing means (thermal sensor) capable 
of sensing the signal and sending it to the signal 
processing system. 

Referring to the drawings and particularly to Figure 

55 1 , the flowmeter system 2 of the present invention gen- 
erally includes a flow transducer 4, a frequency encod- 
ing system 6 and a signal processing system 8. A power 
supply 10 is shown specifically for a heater and is used 
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to provide the necessary energy on instructions f rom the 
frequency encoding system 6 to the flow transducer 4 
to produce a detectable signal (e.g. a thermal pulse) 
. which is convected in the direction of the flow of the fluid 
(e.g. gas). In accordance with the present invention, a 
thermal pulse is generated and supplied to the trans- 
ducer 4 so that it may be carried in the direction of flow 
of the fluid. 

The flow transducer 4 includes a heater 12 and at 
least one sensor, preferably at least one pair of down- 
stream sensors 14, 16. As shown in Figure 1 , the heater 
12 is positioned within the flowpath. It should be noted, 
however, that the heater may be positioned outside of 
the flowpath so long as the heater 1 2 generates a ther- 
mal signal or pulse which is carried downstream in the 
direction of the fluid flow and is sensed in order by the 
downstream sensors 1 4, 1 6. In a preferred aspect of the 
present invention, the transit time or time of flight of the 
thermal pulse between the sensor 14 and the sensor 16 
is used as an indication of a flow characteristic of the 
fluid. 

In accordance with a principle aspect of the present 
inventbn, the frequency encoding system 6 is opera- 
tive ly connected to the heater power supply 1 0 to power 
the generation of thermal pulses which can be sensed 
and processed. In accordance with embodiments of the 
invention discussed hereinafter, the effects of thermal 
diffusion are minimised and the performance character- 
istics of the flowmeter are enhanced. The signals 
sensed by the sensors 14, 16 are processed in the sig- 
nal processing system 8 in a manner which provides 
translation into a desirable flow characteristic consistent 
with the high performance goals described herein. 

The construction of a typical flow transducer em- 
ployed in the present invention can be seen with refer- 
ence to Figures 1 -4. The flow transducer 4 includes a 
housing 20 defining a flowpath 22 which is traversed by 
the heater 1 2 and an optional heater 1 8 for bi-directional 
flows as discussed hereinafter and first and second sen- 
sors 14 and 16. The heater and sensors are suspended 
within the flowpath in desirable locations by securing the 
ends thereof to the housing through the use of a me- 
chanically secure, electrically conductive connection, 
such as a conductive cement. Means of suspending the 
heater and sensors within the flowpath is within the 
knowledge of those skilled in the art. 

The housing of the flow transducer 4 includes an 
inlet 26 and an outlet 28. It will be understood that for 
bi-directional flows, the inlet 26 can, when the direction 
of flow is changed, become the outlet. At the inlet or at 
both ends of a bi-directional transducer there is prefer- 
ably provided devices which serve to condition the flow 
of the fluid as it comes into contact with the heater 12 
and the sensors 14, 16. 

Flow conditioning may be obtained by a variety of 
flow conditioning devices. Such devices can include 
tube bundles and fins. As specifically shown in Figures 
3 and 4, a mesh 30 made of, for example, metal or plas- 


tic is provided to condition the flow of the fluid into and 
out of the transducer 4. Flow conditioning is made nec- 
essary by variations in the geometry of the upstream 
flow conduit (e.g. bends and elbows) which gives rise to 

s distorted flow profiles. Where bends or elbows are 
placed in the path of the fluid flow, there often results a 
distorted flow profile which in turn results in an uncon- 
ditioned or poorly conditioned flow. 

It is desirable to condition the flow of the fluid such 

10 as by the mesh 30 to avoid errors associated with dis- 
torted profiles. Particularly in healthcare applications for 
monitoring patient respiratory flows proximally, velocity 
profile distortion will typically be a problem because of 
the proximity of customary clinical airway fittings such 

is as an elbow and a Y-fitting. 

The mesh 30, preferably positioned substantially 
perpendicular to the longitudinal axis of the conduit, sig- 
nificantly enhances the accuracy of flow measurements 
under disturbed flow conditions by conditioning the flow 

20 prior to the generation of the thermal pulse by the heater 
12. An additional advantage of the mesh 30 is that it 
provides protection for the heaters and sensors con- 
tained within the transducer from accidental contact by 
the user. 

25 The heater 12 and sensors 14 and 1 6 employed in 
the present invention may be any thin element capable 
of radiating thermal energy such as solid wires made 
from metals such as, for example, platinum or nickel. In 
a preferred form of the invention, the heaters and/or sen- 

30 sors are thin film conductors of the type shown in Figure 
5. 

Referring to Figure 5, a film conductor 32 preferred 
for use in the present invention comprises a solid elec- 
trically non-conducting core 34 and a thin metallic film 

35 36 covering at least a portion of the core. A most pre- 
ferred construction of the thin film conductor 32 is an 
optical fibre coated with a thin film of platinum metal. 
The non-conducting optical fibre core 34 preferably has 
a diameter in the range of from about 10 to 40um The 

40 thin metallic film can be applied by vacuum sputtering 
or similar techniques that are capable of providing thin 
metallic films, preferably on the order of from about 500 
to 1000A. 

It is preferred that the electrical resistance of the thin 
45 film conductor be in the range of from about 25 to 50 
ohms/mm. The temperature coefficient of resistance 
(TCR) is a function of the property of the material and 
the thickness of the film. It is preferred that the TCR be 
as high as possible, most preferably at least .001 ohms/ 
50 ohm°C. 

The thin film conductors 32 preferably have a uni- 
form coating along the entire length of the thin film con- 
ductor and a uniform diameter. This construction allows 
thermal pulse generation and sensing performance to 
55 be uniform along the length of the film conductor. It is 
highly desirable that the individual film conductors em- 
ployed in the flowmeter system of the present invention 
have essentially the same or very simitar configuration. 
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This ensures that heater peak temperatures are similar, 
and that sensor thermal responses are the same to pre- 
serve linearity performance and ensure reliability of the 
system for determining the flow direction. 

The signal encoding system 6 as shown in Figure 
1 enables the heater 12 to generate thermal pulses 
which are carried in the direction of the flow of the fluid 
so that they may be sensed by the sensors 14 and 16. 
The shape of the thermal pulse has a significant impact 
on the performance of the time of flight measurement, 
particularly with respect to gas composition independ- 
ence. 

In accordance with the present invention, the heater 
can be completely cooled down prior to the generation 
of the next thermal pulse. Alternatively, particularly 
when the generation of the thermal pulse is through the 
use of a high frequency (e.g. 510 Hz), the heater need 
not cool down to its original starting temperature (e.g. 
ambient) prior to the generation of the next thermal 
pulse. 

What is preferred in accordance with the present 
invention is that the shape of the thermal pulse be sub- 
stantially uniform about the peak temperature thereof. 
As used herein the definition of term "substantially uni- 
form" depends on whether or not the heater completely 
cools down. If it does the term "substantially uniform" 
shall generally mean that the ratio of the response time 
of the heater to reach the peak temperature to the time 
it takes the heater to cool down is from about 1 : 1 0 to 1 0: 
1. The term "substantially uniform" when the heater 
does not completely cool down as it pertains to high fre- 
quency signals shall generally mean that the time it 
takes for the heater to reach the peak temperature is 
from about 20 to 60% of the time between peak temper- 
atures of successive thermal pulses. 

Efforts at measuring the time of flight of a thermal 
pulse are affected by thermal diffusion. The term "ther- 
mal diffusion" is the process by which a temperature dis- 
tribution in a fluid undergoes changes in terms of mag- 
nitude and shape of the temperature distribution over a 
fixed period of time. Thermal diffusion effects are de- 
scribed by the diffusion equation which may be ex- 
pressed in two spatial dimensions as indicated in Equa- 
tions (1)and (2) below. 

2 2 

d{ ax 2 dy 2 

where: 

K= r^ < 2) 

9 is the temperature difference distribution; t is time; x, 
y are spatial co-ordinates; k is the thermal diffusion co- 
efficient of the gas or liquid etc., commonly referred to 


as thermal diffusivity; k is the thermal conductivity; C p 
is the specific heat at constant pressure; and p is the 
density. 

Inspection of the diffusion equation reveals a fun- 
5 damental trend which is that the time rate of change of 
a thermal pulse is more rapid in spatial regions where 
the gradient of temperature change is increasing. Of 
particular importance is the fact that sharply varying 
thermal pulses, and their constituent parts, will tend to 
10 rapidly evolve into smoothly varying thermal pulses 
which will subsequently change shape more slowly. This 
means that the diffusion process not only disperses a 
thermal pulse but it also acts on the shape of the thermal 
pulse and changes it at a rate dependent on the initial 
is shape. 

In order to generate a thermal pulse, it is necessary 
to apply power to the heater to heat it up and then re- 
move the power and allow it to cool down. Prior art sys- 
tems have sought to achieve peak temperatures within 

20 a short time in order to provide clear definition of the 
thermal pulse. This tends to lend a steep and/or straight 
leading edge to the thermal pulse. The initial thermal 
pulse is thus imprinted on the flow of the fluid and during 
the time it takes to arrive at the first sensor becomes 

2S modified due to diffusion. 

The present invention preferably generates thermal 
pulses in which the shape of the thermal pulse is sub- 
stantially uniform about the peak temperature as previ- 
ously described. The signal processing system senses 

30 the entire thermal pulse and does not have to rely on 
leading edge schemes as in the prior art. In the present 
system it is not necessary to wait for the heater to cool 
down between pulses, because the timing is based on 
the shapes of the thermal pulses at the two sensor sig- 

35 nals, and the heater may thus be repetitively pulsed at 
a fixed rate of frequency. A power pulse shape is then 
intentionally selected in order to provide a thermal pulse 
shape which under steady flow conditions is sufficiently 
symmetrical in order that the shape changing effect due 

40 to diffusion produces negligible changes in the shape of 
the thermal pulse. 

In accordance with a preferred aspect of the present 
invention, the heater 12 is pulsed in a manner such that 
complete cool down of the heater does not take place. 

45 in this way, the heater can be more rapidly pulsed, and 
the flow sensor response time can thereby be reduced 
over conventional systems. 

The spatial relationship between the heater and the 
sensors and the sensing of the thermal pulse between 

so the sensors are important considerations in minimising 
the effects of thermal diffusion on time of flight meas- 
urements. As previously indicated the time rate of 
change of a thermal distribution is more rapid where the 
gradient of temperature change over distance is higher. 

55 Thus, a greater change in the shape of the thermal pulse 
can be expected between a heater and a sensor than 
between the more evolved temperature distributions en- 
countered at the spaced apart downstream sensors. It 
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is therefore preferred to determine the time of flight 
measurement between the sensors rather than between 
the heater and the first sensor. 

The position of the heater relative to the sensors 
must also be considered. If the heater is very close to 
the first sensor, the shape of the thermal pulse win un- 
dergo little change as it is carried by the fluid from the 
heater to the first sensor. A more dramatic change in 
shape of the thermal pulse may therefore occur between 
the sensors rather than between the heater and the first 
sensor. It has been determined that the distance be- 
tween the heater and the first sensor is preferably at 
least about 1/3 the distance between the first and sec- 
ond sensors. 

In a preferred form of the invention, the distance be- 
tween the heater and the first sensor is the same as the 
distance between the first and second sensors. 

The time of flight flowmeter measures transit time 
and based on a knowledge of the distance separating 
the sensors, the fluid velocity and not the volumetric 
flowrate is initially measured. The conversion of fluid ve- 
locity to volumetric flowrate can lead to errors. It is 
known in fluid dynamics that for any particular set of flow 
conditbns, the velocity within a flow conduit varies ac- 
cording to position within the conduit. As a result, the 
centre position does not normally reflect the mean value 
or provide a straight line response. Therefore, position- 
ing the sensors in the centre position (i.e. on or across 
the longitudinal axis of the flow conduit) of the conduit 
can lead to errors. 

In addition, the flow profile of the fluid passing 
through a conduit is known to be a function of the flow- 
rate and the composition of the fluid. The most dramatic 
variations in flow profile accompany a change from lam- 
inar flow to turbulent flow Thus, when the flow of a fluid 
undergoes a change from laminar flow to turbulent flow, 
significant errors in the determination of volumetric flow- 
rate can result. 

In accordance with the present invention, both lam- 
inar and turbulent flows can be measured in the same 
flowmeter to obtain accurate volumetric flowrates. It has 
been determined that the employment of a line sensor 
to sense a detectable signal (e.g. a thermal pulse) along 
a line traversing the direction of flow of the fluid reduces 
errors due to changes in flow regimes, especially when 
the line sensor is off -axis (i.e. is not on or does not cross 
the longitudinal axis of the conduit). 

An example of a line sensor and its off-axis position 
in the conduit is shown in Figures 3 and 4. The line sen- 
sor is capable of sensing a thermal pulse along at least 
a substantial length of the sensor. Line sensors are dis- 
tinguished from point sensors which sense thermal puls- 
es at a particular point in the direction of flow of the fluid. 
It will be noted that the sensors 1 4 and 1 6 are positioned 
off the longitudinal axis of the conduit (i.e. they are not 
on and do not cross the longitudinal axis) in a line tra- 
versing the flow direction of the fluid. The preferred off- 
axis location depends on the shape of the conduit. For 


example, for a conduit having a circular cross-section, 
the location of a line sensor is determined by the Equa- 
tion (3) 

r/R = 0.4 (3) 

Wherein r is the radial distance of a line sensor from the 
centre of the conduit and R is the radius of the conduit. 
10 The selection of a frequency for the generation of 
thermal pulses is an important criteria for achieving high 
performance of the flowmeter. The ability of the signal 
processing system to discriminate the thermal pulses in 
order to make a time of flight measurement, is related 
is to the relative phase angle of the signal. For example, 
a time of flight of a half cycle represents a relative phase 
angle between the signals at the first and second sen- 
sors of 180°. This is a relatively easy measurement 
compared with a smaller relative phase angle such as 
20 1°. Thus, the ability to resolve differences in the time of 
flight, which directly impacts on the accuracy of the 
measurement, is fundamentally affected by the choice 
of operating frequency. In a preferred form of the present 
invention, the signal encoding system generates a sig- 
25 nal having the combination of simultaneous multiple fre- 
quencies including at least one higher frequency and at 
least one lower frequency. Higher frequencies are pre- 
ferred when a faster response time is desired. Lower 
frequencies are preferred because they can more read- 
so ily be used to measure the time of flight of thermal pulses 
at lower flowrates. 

In accordance with the present invention, the em- 
ployment of a multiple frequency signal can extend the 
performance level of the flowmeter as compared with a 
35 signal frequency system. The employment of a bwer 
frequency allows the system to operate at lower flow ve- 
locities while higher frequencies achieve greater accu- 
racy through improved repeatability as well as faster re- 
sponse times. 

40 The selection of desirable higher and lower fre- 
quencies is made with regard to the range of flowrates, 
dimensions of the transducer, and thermal process lim- 
itations encountered in the flowmeter. By way of exam- 
ple and especially for gases, typical higher frequencies 

45 can be in the range of from about 200 to 700 Hz (e.g. 
510 Hz) while typical lower frequencies are in the range 
of from about 10 to 50 Hz (e.g. 33 Hz). In the case of 
liquids the operating frequencies normally tend to be 
lower than for gases. In operation, the signal comprised 

50 of the higher and lower frequencies is employed to gen- 
erate a thermal pulse. The sensing of that signal and its 
processing is made by selectively filtering to obtain the 
higher and lower frequency components. Details of the 
processing of the multi-frequency signal will be dis- 

55 cussed hereinafter. 

In a preferred form of the invention the thermal 
pulse rate is varied at a high frequency which is frequen- 
cy modulated or frequency and amplitude modulated. If 
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the heater voltage is pulsed at a rate which is frequency 
modulated, the power spectrum of the fluid convected 
thermal pulse can evidence the FM spectrum, which in- 
corporates the carrier frequency and FM sidebands plus 
the lower, modulation frequency. This is because the flu- 
id-heater combination partially demodulates the fre- 
quency modulated power wave form at the heater as a 
result of thermal inertia, and both frequencies are trans- 
ferred to the thermal pulse and convected simultane- 
ously. An example of such a system is shown in Figure 6. 

Referring to Figure 6, there is shown an arrange- 
ment of a system by which a suitable frequency and am- 
plitude modulated thermal pulse is generated. An FM 
wave form generator 40 generates a thermal pulse, pref- 
erably in the shape of a square with alternating polarity 
shown by numeral 41 . This wave form is then continu- 
ously fed into a power supply 42 which provides the nec- 
essary voltage for the heater to provide thermal pulses 
at the desired frequency. At the same time or alternately, 
a trigger pulse is generated once every FM cycle and 
sent to an AM function generator 44. The resulting am- 
plitude modulated wave form shown by numeral 43 is 
sent to the power supply 42. The voltage output level is 
varied in alternating increments once every FM cycle 
period e.g. 30 milliseconds is the FM frequency is 33 
Hz, which thereby generates a step shaped amplitude 
modulation to the heater (e.g. 16.5 Hz). 

An example of a signal power spectrum generated 
by an FM encoded heater is shown in Figure 7. Referring 
to Figure 7, there is shown a power spectrum of the up- 
stream sensor signal for the gas convected thermal 
pulse when the centre frequency is 510 Hz modulated 
at 33 Hz with an index of modulation of 1.5. The gas 
convected thermal pulse contains the FM spectrum, 
which in this example includes the centre frequency and 
two sidebands. 

It is possible to generate additional frequencies us- 
ing dual tone FM modulation. This can be done with 
mixed modulation frequencies of, for example, 33 Hz 
and 10 Hz. 

As discussed hereinafter, control of the peak tem- 
perature of the heater is advantageous and is achieved 
by controlling the voltage and hence power to the heater. 
Due to variations in heater cooling rates, as the gas ve- 
locity or composition varies, the peak temperature of the 
heater varies if the power to it is not controlled. Ampli- 
tude modulation is implemented as part of the system 
for controlling the peak temperature of the heater. Con- 
trol of the power to the heater is therefore integrated with 
the frequency encoding system. 

Employing an FM encoded power wave form, ad- 
justments to the heater power cannot be made at the 
time of every high pulse rate frequency (e.g. 510 Hz) 
without impairing transfer of the modulation frequency 
(e.g. 33 Hz). Therefore, adjustments to the pulse voltage 
and hence heater power are made once every FM cycle 
(e.g. 30 milliseconds) as in fact indicated by the trigger 
pulse to the AM generator as shown in Figure 6. At that 


time small errors in the adjustment will generally give 
rise to an amplitude modulation at half the FM modula- 
tion frequency (e.g. 16.5 Hz). Therefore, the preferred 
heater voltage wave form is of the type wherein the AM 
s component is at the lowest frequency such that the AM 
frequency is equal to one half the FM frequency. All 
three voltage wave forms (FM, FM + AM and AM) are 
alt viable heater voltage wave forms for real-time oper- 
ation without phase ambiguity. 
10 The resulting AM modulation is preferably control- 
ted by applying known alternating voltage differences at 
each voltage adjustment. In this way a square wave am- 
plitude modulation is reliably generated at a known 
depth of modulation. In another system amplitude mod- 
's ulation can also be used to generate the lower frequency 
instead of frequency modulation and this alternative is 
particularly important if the fluid being measured is a liq- 
uid. In this case the liquid will carry the heat from the 
heaters so rapidly that FM demodulation would not oc- 
20 cur and amplitude modulation would have to be used. 
A number of viable options therefore exist for gen- 
erating the lower frequency components of the simulta- 
neous multifrequency signal. These include FM or AM, 
or FM and AM. 

25 Having generated a detectable signal (e.g. a ther- 
mal pulse), preferably having multiple frequency com- 
ponents, a signal processing system is provided to 
sense the detectable signal and to provide a flowrate 
measurement meeting the performance criteria as de- 

30 scribed herein. 

A preferred signal processing system is shown in 
Figure 8. A flow transducer 4 contains sensors (not 
shown) as previously described in connection with Fig- 
ure 1. The sensors are connected via a pair of shielded 

35 leads (channels 1 and 2) to the inputs of a pair of AC 
coupled high gain amplifiers shown by numeral 50. 

The sensors are resistive devices which change re- 
sistance approximately linearly with temperature. 
Therefore, the amplifier 50 is required to supply a stable 

40 DC excitation current. As previously indicated, the sen- 
sors are preferably thin film conductors having a non- 
conducting core and a thin film of metal as shown in Fig- 
ure 5. As the resistance of the thin film conductors 
changes in response to the presence of a thermal pulse, 

45 a signal voltage is developed. The change in resistance 
with temperature is dictated by the temperature coeffi- 
cient of resistance of the sensor which is equal to or re- 
lated to the temperature coefficient of resistance of the 
sensor film material. 

50 The excitation currents for the two sensors are pref- 
erably within 5% of each other. Within this range, the 
signal voltages at the two sensors, which are propor- 
tional to the thermal pulse magnitude may be reliably 
compared for the purpose of determining flow direction. 

55 By way of example, the preferred excitation currents us- 
ing the preferred platinum thin film conductors previous- 
ly described is in the region of 2 mA. The preferred ex- 
citation currents are in the range of from about 0.4 to 1 2 
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mA with a most preferred range of from about 1 to 4 mA. 

After the front end amplifiers 50, the strength of the 
channel 1 and channel 2 signals are preferably in the 
range of a few tenths of a volt peak to a peak. Low fre- 
quency components tend to have a higher magnitude 5 
than the higher frequency components. In order to utilise 
common signal pathways, it is preferred to even out the 
magnitudes. The signals are coarsely partitioned into 
low and high ranges and sent to respective tow and high 
frequency branches. In the preferred system, the chan- 10 
nel 1 and channel 2 tow frequency signals tend to be up 
to 6 times greater in magnitude than the high frequency 
signal and are sent to an attenuator 52 to reduce their 
magnitude relative to the high frequency signals by a 
factor of 6:1. The pairs of signals in both the tow and is 
high frequency branches are sent through respective 
partitioning filters 54 which are preferably digital filters 
such as switched capacitor devices. Cut off frequencies 
must be selected such that at the signal frequencies of 
interest, group delay differences are very small com- 20 
pared with the time-of-flight measurement to thereby 
avoid systematic timing errors. 

The resulting signals are sent to respective non-lin- 
ear amplifiers which amplify the signals by identical fac- 
tors over most of the voltage input range to provide a 
range matching to the analogue to digital converters 58. 
Prior to entering the analogue to digital converters, the 
low and high frequency signals are recombined in a pair 
of summing amplifiers 60 and then converted into digital 
format via the analogue to digital converters 58 for fur- 30 
ther signal processing via the digital signal processing 
device 62. 

The flow of data into the digital signal processing 
device 62 may be at a rate up to 64 kHz but is preferably 
1 6 kHz. When entering the signal processing device 62, 35 
the signals are first filtered to select specific spectral re- 
gion pairs for subsequent cross correlation as shown 
specifically in Figure 9A. 

Referring to Figure 9A, the channel 1 and channel 
2 signals are spectrally filtered by three pairs of filters in *o 
order to select three specific signal frequency regions 
which span the specific frequencies generated by the 
frequency encoding system for cross correlation over 
the corresponding specific time periods. The first pair of 
filters are, for example, 8 pole, low pass Butterworth 45 
types with a -3 dB 70 Hz cut off frequency. The second 
pair of filters are, for example, 8 pole elliptical types with 
-3 dB points at 30 Hz and 180 Hz and a group delay 
which varies from about 40 ms at 33 Hz to around 12 
ms at 66 Hz. The third pair of filters are 8 pole Cheby- so 
shev types with -3 dB points at 478 Hz and 550 Hz and 
a group delay which varies from about 10 ms at 510 Hz 
to about 20 ms at 477 and 543 Hz. This pass band thus 
encompasses the 510 Hz FM spectrum. 

Group delay is a characteristic of filters and impacts ss 
important timing factors in the transit time measure- 
ment. The typical group delay and signal frequencies 
determines the time at which the channel 1 and 2 signals 


are available for subsequent analysis and time-of-flight 
measurement. It does not determine the minimum time- 
of-flight that can be measured but rather the lag time 
between a change in flow and a corresponding electron- 
ic measurement. This relates to the speed at which the 
filters can pass signals and hence the response time of 
the filters which is indicative of the overall response of 
the system to changes in flow. 

It is a characteristic of the processes by which ther- 
mal pulses are impacted and then conveyed by the gas 
flow that lower frequencies are more effective for signal 
generation at lower flow velocities than higher frequen- 
cies. Hence the use of lower frequencies extends the 
operating range to lower flowrates where higher fre- 
quency signals tend to decay and disappear 

As shown in Figure 9A, the channel 1 and channel 
2 analogue inputs are transferred to the analogue to dig- 
ital converter 58 (see Figure 8) and are thereafterfiltered 
by one of three signal filter pairs. The number 1 or low 
signal filter and the number 3 or high signal filter are 
continuously analysed to determine the peak signal 
strength within the signal period. Hence for low signals 
the peak value is determined within a given set time pe- 
riod (e.g. 60 ms) and the high signals the peak values 
are determined during a second preselected period (e. 
g. 1 .96 ms). Truth values are then generated whereby 
if the signal exceeds certain preselected peak values (e. 
g. 20 dBV below normal operating values), they are con- 
sidered to be viable and valid thermal pulse signals. 
Hence the outputs of these measurements consist of 
four flags with assigned logical truth values representing 
low channel 1 valid (1 ,0); low channel 2 valid (1 ,0); high 
channel 1 valid (1,0); and high channel 2 valid (1,0). 

The three pairs of filter signal data streams feed into 
three simultaneous cross correlation routines called 
LOW, MID and HIGH. The three routines operate on 
synchronous data. The synchronous lines connect LOW 
and HIGH cross correlations with the LOW and HIGH 
signal strength measurements. The signal strengths are 
therefore synchronously measured as the cross corre- 
lations are performed. The MID routine will therefore be 
automatically in phase with the LOW routine. The syn- 
chronisation is in order that the LOW valid and HIGH 
valid flags, always correspond to the same data operat- 
ed on in the current cross correlation operations. 

The LOW routine cross correlates the signals from 
the pair of low pass filters and provides time-of-flight in- 
formation designed for the low flow region of the flow 
dynamic range. For example, if the LOW routine oper- 
ates at 2 kHz, and the analogue to digital data stream 
is operated at 1 6 kHz, the signal prefilter works on every 
eighth data point. The cross correlation base period in 
this example is 60 ms, which corresponds to the period 
of the lowest frequency 16.5 Hz signal component. 

Referring to Figure 9B there is shown a graph show- 
ing a cross correlation function between two sensor sig- 
nals where the signals contain the two frequencies 510 
Hz and 33 Hz. The peak of the function which indicates 
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the time delay between the signals is shown at 2.8ms. 
This time delay is greater than the cycle period of the 
higher frequency (510 Hz) component and thereby dem- 
onstrates the ability of this technique to avoid phase am- 
biguity. 5 

The MID routine can operate, for example, at 4 kHz 
but preferably at 8 kHz. The outputs of the two pairs of 
filters (signal filter pair number 2 and number 3) are add- 
ed together and submitted to the MID cross correlation 
routine which correlates the combined signals of the cy- 10 
cle period of the lowest frequency component (33 Hz) 
which is 30 ms. The HIGH routine cross correlates the 
signals from the pair of band pass filters number 3 and 
is operated at 16 kHz. The base time period is about 2 
ms corresponding to the 510 Hz signal frequency. 

Figure 1 0 shows the signal strength process algo- 
rithm employed in Figure 9A for providing flow zero, flow 
direction, extreme low flow and signal valid outputs. As 
previously indicated, the outputs shown in the digital da- 
ta flow of Figure 9A consist of four flags with assigned 
logical truth values. Figure 10 depicts the subsequent 
signal strength process algorithm which continually op- 
erates on these values to determine flow direction and 
the like. 

The flow zero flag is set to 1 if the signal strength in 
the lowest frequency region is below a predetermined 
signal valid level at both sensors. When the flow zero 
flag is set, the system is designed to respond with a flow- 
rate that equals zero output. Since zero cannot be 
measured, this means instrument zero flow to within the 
low flow sensitivity of the device. The low valid flag is 
assigned a value of 1 when the signals in both channels 
in the lowest frequency region are valid and signifies that 
low flow time-of -flight flow measurements can be report- 
ed. 

The HIGH valid flag is set at 1 such that the signals 
in both channels in the highest frequency region are 
above a predetermined level and that the time-of -flight 
measurements based on high can be used. The forward 
flag is employed for bi-directional operation when two 
heaters are employed. This flag is set to 1 if the signal 
strength at sensor 1 is valid and greater than the signal 
strength at sensor 2. The flag is set to zero if the signal 
strength at sensor 2 is valid and greater than the signal 
strength at sensor 1 . Flow direction is hence determined 
even at low flows when the signal strength of the down- 
stream sensor is insufficient for an appropriate time-of- 
flight measurement at the designated flowrate. 

An optional "extreme low* flag is set at 1 when the 
flowrate is so low that the downstream sensor signal 
strength is insufficient for a valid time-of-flight determi- 
nation but the upstream sensor signal is still valid. The 
flowrate under these conditions is therefore between ze- 
ro and the time-of-flight measurement sensitivity. An es- 
timated value midway between these values may be se- 
lected for output and will provide better flow information 
than a zero output although the accuracy is reduced 
compared with normal operation. This option to the 


present invention provides a greatly enhanced low flow 
sensitivity at reduced accuracy and is particularly desir- 
able in certain clinical neonatal situations. 

Referring again to Figure 9 A, the time-of-flight de- 
termination has available the current flow direction from 
the signal strength routine. The cross correlation rou- 
tines and subsequent peak finds determine the relative 
shift between the two signals over the corresponding 
base period of the signal frequency. If the initially as- 
sumed flow direction, as characterised by the selectbn 
of S A and S 2 data streams, was in the opposite direction, 
and the time-of-flight is corrected at this point by sub- 
tracting it from the base time period, then the channel 1 
and channel 2 designations for all the cross correlation 
routines are switched. 

The individual time-of-flight results of the continu- 
ous cross correlation routines are submitted to data fil- 
ters. The time-of-flight values are first filtered and the 
results are then converted into equivalent flowrates. 
Two types of digital data filters may be employed. The 
first type is called a median filter which operates on the 
previous N data points where N is an odd number and 
is described as having a rank R where N = 2R + 1 . This 
type of filter completely eliminates spurious, isolated da- 
ta points which lie beyond the normal random spread of 
the current data. 

The second type of filter is often referred to as a box 
car filter. This type of filter simply averages the previous 
N data points and hence imposes a 0-100% response 
time on the output signal equal to the time taken for the 
N data points to pass completely through the filter. 

The low routine time-of-flight output is preferably fil- 
tered using only a rank 1 median filter. The MID routine 
is preferably filtered using a rank 2 median filter ahead 
of a 2 point box car filter. The HIGH routine is preferably 
filtered using a rank 7 median filter ahead of a 25 point 
box car filter. 

A slower speed of response but simpler alternative 
to the above system is to utilise only a single cross-cor- 
relation routine which incorporates two or more frequen- 
cy components of the signal in the same manner as the 
MID routine such that the correlation period is the same 
as the period of the lower frequency component. This 
routine has the advantage that it provides a seamless 
technique for measuring transit times that traverse the 
region of phase ambiguity associated with the higher 
frequency signal component. 

Another method for measuring the time-of-flight 
employs a zero crossing timer as shown specifically in 
Figure 11. Here the signal inputs are AC coupled and 
preferably sinusoidal. Simple comparator circuits are 
used to detect the point at which the first signal crosses 
the zero volt level and this event triggers a time counter 
device. The timer is then stopped and the number of 
counts read when a second zero crossing occurs, this 
time on channel 2. 

Referring again to Figure 9A, a main flow algorithm 
is employed to convert all of the data in the digital flow 
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sequence into a flowrate and flow direction output. An 
example of the main flow algorithm is shown in Figure 
12. 

The algorithm initialises and checks thai the signal 
strength and cross correlation routines are running. 
Then, depending on flag truth values from the signal 
strength routine indicating signal validity, the main algo- 
rithm proceeds to select the appropriate flowrate output 
value via a continuously and rapidly looping signal va- 
lidity and result availability check and a transition level 
scheme. 

The basic requirement is to obtain a selection of the 
appropriate correlation routine flow values as the flow- 
rate varies. Hence, if the flowrate according to the low 
routine exceeds a certain lower level then it is more ap- 
propriate to use the result from the MID routine if it is 
available. However, if the MID routine value is greater 
than a certain threshold value flowrate, then it may be 
more appropriate to select the output of the high cross 
correlation routine if it is available. If the flowmeter is 
operational but there is no flow through the transducer, 
then a very fast loop will continually check this condition 
and report a flow of zero. If the flow is instead above 
zero, but nevertheless insufficient to render the low sig- 
nals valid at both sensors, then another rapid loop will 
report the extreme low flowrate estimate and flow direc- 
tion. This provides a fast response, extra low flow sen- 
sitivity which is particularly useful in rapid breathing ne- 
onates in a neonatal application. 

When the low signals are valid, the result from the 
low cross correlation flow measurement is always avail- 
able. This assures that there will be no time delay in the 
initial output as the flow accelerates from zero. Depend- 
ing then on the outputs of the MID and high cross cor- 
relations vs. the selected threshold value flowrates, the 
flowrate value may or may not be subsequently reas- 
signed in traversing the remainder of the algorithm to 
the update output stage at which point control loops 
back to the beginning. 

The heaters and sensors employed in a flowmeter 
are physically close to one another. Because the heater 
is subjected to repeated and relatively rapid applications 
of power, a signal appears at the sensor at the same 
time as the power is applied. This signal, referred to as 
"crosstalk interference* must be distinguished from the 
thermal pulse signal in order to obtain accurate flow 
measurements. Crosstalk interference can lead to po- 
tential timing errors and also adversely affect the sensor 
signal amplifiers. 

The crosstalk signal is generated due to the capa- 
citative and inductive coupling which unavoidably oc- 
curs between unshielded conductors such as heaters 
and sensors. 

The magnitude of the interference signal is deter- 
mined by the geometry of the transducer, the rate of in- 
terfering voltage or current application at the heater and 
the sensitivity of the sensor as compared with the ther- 
mal signal sensitivity. The magnitude of the interfering 


signal is very sensitive to the steepness of the power 
pulse. The magnitude of this interfering signal can be 
alleviated somewhat by smoothing filters applied to the 
heater voltage waveform. In the present system, if the 

s crosstalk interference signal were incorporated into both 
sensor signals, there would be distortion and possible 
timing errors. 

In accordance with the present invention, this prob- 
lem can be eliminated by alternating the polarity of suc- 

10 cessive power pulses. The fluid, which has no appreci- 
ation for the direction of current flow in the heater, re- 
ceives the thermal pulse at the pulse rate frequency, 
whereas crosstalk interference is created at exactly half 
the pulse rate frequency Thus, there is no crosstalk in- 

15 terfering spectral component at the pulse rate frequen- 
cy. Accordingly, the sensor signals which are subse- 
quently electronically filtered, are eventually completely 
devoid of crosstalk interference and no timing errors oc- 
cur. 

20 in fluid flow sensing applications the thermal time- 
of-flight peak heater temperatures tend to be relatively 
high. In some applications, particularly those involved 
in healthcare, federal regulations provide strict limits on 
the maximum operating temperature of healthcare ap- 

25 paratus. The present invention therefore provides a 
heater temperature control system which works in con- 
junction with the frequency encoding system previously 
described in connection with Figure 1 . 

During normal operation of the thermal time-of-fiight 

30 transducer, as the flowrate varies, the rate at which heat 
is lost by the heater also varies. The rate of heat transfer 
from the heaters to the fluid flow depends not only on 
the velocity or rate of flow, but also on the composition 
of the fluid. If the heater pulse power is fixed at a preset 

35 level, then flowrate and composition variations will pro- 
duce heater peak temperature variations. Since the dif- 
ference between heat loss rates at the extremes of the 
range of flows through the transducer can be substan- 
tial, there are likewise substantial variations in peak 

40 heater temperatures. The heater temperature control 
system of the present invention automatically and ac- 
tively controls the peak heater temperature to prevent 
heater material failure, ignition of any flammable gases 
and performance variations. 

45 it should be also noted that in the healthcare indus- 
try, especially in a critical care setting to measure respi- 
ratory flows of patients under anaesthesia, federal reg- 
ulations require that the surface temperature of compo- 
nents in the airway not exceed 300°C. Therefore, the 

so heater temperature control system of the present inven- 
tion may be designed to actively control the peak tem- 
perature of the heater so that it does not exceed 300°C. 
The present invention may be also provided with the 
means for controlling the peak temperature to a set point 

55 (e.g. between 100°C and 200° C) suitable for providing 
practical signal strengths and good performance. By 
employing the heater temperature control system of the 
present invention, preferred FM signal encoding is pre- 
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served and the heater power waveform may be ampli- 
tude modulated at half the FM frequency. 

The heater control system varies the heater voltage 
such that during normal operation the strength and 
spectral integrity of the thermal pulse signals are pre- 5 
served. The system generally comprises two linked tem- 
perature control subsystems which utilise peak heater 
temperature measurement information obtained via an 
electronic heater resistance measurement. The first 
subsystem is a limiter circuit which is set to continuously io 
and automatically ensure that the peak temperature (e. 
g. 300°C) is never exceeded. It operates by removing 
the power applied to the heater during a thermal pulse 
whenever the heater resistance exceeds a preset value. 
The first subsystem therefore overrides the second sub- '5 
system in the event that the peak heater temperature is 
exceeded. 

The second subsystem comprises a feedback con- 
trol loop in which a set point temperature is pre-pro- 
grammed which is less than the pre-set maximum tern- 20 
perature defined above for the first subsystem. This set 
point temperature represents desired normal operating 
heater peak temperature. 

The peak heater temperature is measured, com- 
pared with the set point temperature and the heater volt- 
age adjusted in accordance with a control algorithm 
once every FM cycle period. In addition, the heater volt- 
age adjustment is further used to provide a lower fre- 
quency component to the signal via amplitude modula- 
tion as part of the signal encoding system. 

Referring to Figures 13 - 15, the heater control sys- 
tem 80 includes a voltage divider circuit 82 to measure 
the heater temperature by measuring the resistance of 
the heater (e.g. thin film conductor). The voltage divider 
circuit is more specifically set forth in Figure 14. Refer- 
ring to Figure 14, the heater represented as a resistive 
element R H forms part of a potential divider with another 
fixed value resistor, R B , with the power supply connect- 
ed thereto. The potential across the fixed base resistor 
is then amplified. It is therefore possible by knowing the 
input voltage, Vj, the divider voltage, V B , the gain, G, 
and the base resistance value to routinely determine the 
resistance of the heater. For example, the R H for a single 
heater which can be employed in the present invention 
may be in the region of about 400 ohms. The base re- 
sistor has a smaller value and is a stable, adequately 
power rated resistor of about 5 ohms while the gain is 
about 10. 

Referring again to Figure 13, the divider voltage is 
connected to a sampling circuit 84. The sampling circuit 
samples the value of the divider voltage at specified 
times preferably corresponding to the heater peak tem- 
perature and the FM modulation frequency, cycle peri- 
od. 

The sampling circuit 84 has three inputs and gen- 
erates two outputs. As described in connection with Fig- 
ures 13 and 15, the circuit accepts the FM modulation 
waveform 85 from the FM function generator 86. The 


cycle period (e.g. 33 Hz) modulation frequency is de- 
tected by the positive (as shown specifically in Figure 
15) or negative zero crossing. When this condition is 
flagged the next power pulse is detected via the power 
waveform input 95 and the leading edge starts a timer 
preset to time out at a time just less than the width of 
the pulse as shown in Figure 1 5. Clearly the temperature 
of the heater will tend to reach a maximum at a time 
corresponding to the trailing edge of the power pulse. It 
is preferred to sample the temperature just prior to that 
point. 

At this time the output 87 of the voltage divider cir- 
cuit 82 is sampled in the sampling circuit 84. The sam- 
pled voltage 91 is sent as an input to the heater temper- 
ature control algorithm 89. Simultaneously a synchro- 
nising pulse 93 is generated which has a width substan- 
tially less than the FM cycle period and is about, for ex- 
ample, 10 ms wide and is also sent to the A/D and D/A 
devices which mediate the flow of data in and out of the 
control algorithm 89. The trailing edge of the synchro- 
nising pulses is used to strobe the A/D and D/A devices. 
Thus the sampled voltage divider voltage is read once 
by the control algorithm as a new value each time the 
A/D device is so activated. 

An example of a temperature control algorithm for 
use in the heater control system is shown in Figure 16. 
When the limiter circuit (see Figure 1 3) is active the out- 
put V variable is adjusted to a predetermined non-zero 
base value. This base value actually corresponds to a 
low pulse voltage setting, preferably in the range of from 
about 5 to 10 volts. At this level, the heater peak tem- 
perature is close to the set point peak temperature at 
zero flow. 

The flow chart shows that there is always control 
and knowledge of the pulse voltage setting. Given there- 
fore the sampled divider voltage and the predetermined 
values for the base resistance of the divider, the gain of 
the divider circuit, the ambient or cold temperature value 
of the heater resistance and the TCR value of the heater, 
the heater temperature is calculated and the necessary 
voltage adjustment to bring the peak heater temperature 
closer to the targeted temperature is determined. The 
new output voltage is written to the location RANGE 
CONTROL. When the output device (e.g. a digital to an- 
alogue converter) is strobed by the negative edge of the 
SYNC PULSE this value is converted to an appropriate 
analogue voltage via the digital to analogue converter 
which is connected to the SET VOLTAGE input of the 
power supply. This input voltage linearly controls the 
pulse voltage via a range factor which is built into the 
power supply such that the voltage output range of the 
digital to analogue converter ratios the output voltage 
range of the power supply. 

The feedback temperature control system adjusts 
the heater power pulse voltage once every FM cycle pe- 
riod. This timing is necessary in order not to destroy the 
FM demodulation thermal pulse signal component. 
Therefore, the control algorithm adjusts the pulse volt- 
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age every FM cycle period. At this point, amplitude mod- 
ulation is also generated. The AM modulation may be 
applied prior to writing the new output voltage variable 
to the RANGE CONTROL LOCATION. This can be ac- 
complished by applying alternately positive and nega- 5 
tive voltage increments to the voltage output value. 

Other amplitude modulation waveforms and fre- 
quencies may be generated at this point. Amplitude 
modulation can be solely used to generate the lower fre- 
quency components) of the signal for fluids in general, 10 
especially when the fluid is a liquid. This can be accom- 
plished in the simple step waveform manner previously 
described. If instead a more complicated multifrequency 
amplitude modulation is required, a continuous ampli- 
tude modulation function could be applied by an external is 
amplitude waveform generator. The heater peak tem- 
perature control could then be achieved by instead ap- 
plying incremental adjustment to the mean of that func- 
tion. 

Temperature control as discussed previously can 20 
be applied to a bipolar, bi-directional configuration em- 
ploying two pairs of heaters and sensors so that four 
heaters are present in the transducer The peak temper- 
ature of multiple heaters in a single transducer can be 
controlled when the heaters are essentially the same 25 
and are electrically connected together. The overall am- 
bient resistance values for the combination are then de- 
termined and entered into the control program in the 
same way as previously described along with the TCR 
value. The heater resistances should then be similar in so 
order that the individual peak temperatures are suffi- 
ciently close to each other. 

The transducer employed in the present invention 
may optionally include multiple combinations of parallel 
pairs of heaters and sensors, preferably two pairs locat- 35 
ed directly opposite one another, preferably equally dis- 
tant from the axis at the preferred off -axis position. This 
multiplanar flowmeter can provide increased accuracy 
by enhancing repeatability due to increased signal 
strength and decreased noise. Such a system is also 40 
more highly insensitive to potential flow profile asymme- 
try errors by reducing the demand for flow conditioning. 

The flowmeter in this embodiment comprises a plu- 
rality of signal generating and sensing units, with each 
unit comprising at least one signal generator for passing *5 
the detectable signal into the fluid and at least one signal 
sensor for sensing the signal downstream of the point 
when the detectable signal is passed into the fluid. Each 
of the units is positioned in a plane along the mean flow 
direction of the fluid, said planes being parallel to each so 
other so that the units are spaced apart and parallel to 
each other in the mean flow direction of the fluid. 

The flowmeter may be provided with more than one 
combination of heaters and sensors. This biplanar ar- 
rangement provides improved accuracy by more clearly ss 
differentiating between the desired signal and undesired 
noise. The biplanar arrangement is also less sensitive 
to errors arising from flow profiles and therefore requires 


less flow conditioning to obtain accurate flow measure- 
ments. 

Referring to Figure 17, there is shown a first com- 
bination 90 of heaters and sensors including a heater 
92 and respective sensors 94, 96. There is also provided 
a second combination 98 including a heater 100 and re- 
spective sensors 102 and 104. The respective combi- 
nations 90 and 98 of heaters and sensors are positioned 
directly across from each other (i.e. in the same longi- 
tudinal position within the flow path). 

The adjacent sensor pairs 94, 96 and 102, 104 are 
preferably serially connected so that excitation currents 
are exactly the same and the heaters 92, are all oper- 
ated from the common heater power supply, so that ther- 
mal pulse generation in the two planes is simultaneous. 
Assuming, for example, equal heater and sensor resist- 
ances, and hence sensitivity, the signal magnitude at the 
sensor pairs is now doubled. However, if the noise com- 
ponent of the signals is randomly generated, then the 
noise instead only increases by 72 . Hence, the signal 
to noise ratio, which along with other factors determines 
the repeatability performance, should increase by about 
41%, and 

_2 
72 

and has been found to increase by as much as 100%, 
which provides enhanced repeatability performance. 
The high rate of improvement indicates that the intensi- 
ties of some components of the velocity fluctuations 
which give rise to much of the 'noise' are themselves 
anticorrelated and cancel. 

With the use of a multiplanar arrangement, the de- 
tectable signals can also be autocorrelated at adjacent 
sensor pairs and then the autocorrelations subsequent- 
ly crosscorrelated. The time of flight determination then 
involves cross correlating the pair of autocorrelation 
functions. This results in an even greater enhancement 
in the signal to noise ratio and hence superior accuracy 
through enhanced repeatability at a given speed of re- 
sponse throughout the flow range. 

The preferred signal processing system for use with 
a multiplanar transducer is shown in Figure 18. In this 
arrangement, the signal processing system connects 
the sensors individually to four separate channels des- 
ignated 1 A, 1 B, 2A and 2B. Signals from the four chan- 
nels are sent to the front end amplifiers 70, the signal 
partitioning 72 and the signal filters 74 operate as de- 
scribed in connection with the signal processing system 
of Figure 7. The components 70, 72 and 74 are dupli- 
cated to provide capacity for extra channels. 

The signal strength measurements can be modified 
to include the A and B signals as the first algorithmic 
step to obtain average channel 1 and channel 2 quan- 
tities whereupon the subsequent process remains es- 
sentially the same as that described above for Figure 8. 
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Autocorrelation as previously described is essentially 
the same as cross correlation with the difference that 
the two signals entering the correlation are considered 
different representations of the same signal. 

The channel 1 and channel 2 autocorrelation func- 
tions are generated between each pair of plane A and 
B sensor signals for each channel and for each signal 
frequency range as defined by the signal filters 1 , 2 and 
3. The autocorrelation functions are then submitted to 
the cross correlation function routines to be processed 
in the same way as if they were channel 1 and channel 
2 filter signals as in the system of Figure 8. It will be 
noted, however, that the actual numbers will be much 
larger and could therefore require factoring to avoid 
overflow in the signal processing. 


Claims 

1 . A device for measuring a flow characteristic of a flu- 
id comprising: 

(a) at least one signal generating means 12 po- 
sitioned in proximity to the flowpath of the fluid 
for generating a detectable signal and for deliv- 
ering the detectable signal within said flowpath; 
and 

(b) at least one pair of longitudinally spaced 
apart sensors 1 4, 1 6 positioned downstream of 
the signal generating means 1 2 for sensing the 
detectable signal as a measure of the time of 
flight of the detectable signal between the pair 
of sensors. 

2. A device tor measuring a flow characteristic of a flu- 
id comprising: 

(a) at least one signal generating means 1 2 po- 
sitioned in proximity to the flowpath of the fluid 
for generating detectable signals and for deliv- 
ering the detectable signals within said flow- 
path; and 

(b) at least one sensor positioned downstream 
of the signal generating means 12 for sensing 
the detectable signals as a measure of the time 
of flight of the detectable signals. 

3. A device as claimed in claim 1 or 2, in which the 
signal generating means 12 generates detectable 
signals each in the form of a thermal pulse having 
a uniform shape about a peak temperature. 

4. A device as claimed in claim 3, in which the signal 
generating means is a heating means 12 including 
a source of power 10 for generating a thermal pulse 
and means for controlling the generation of thermal 


pulses to prevent the heating means 12 from ex- 
ceeding a preselected maximum temperature. 

5. A device as claimed in claim 4, in which means is 
5 provided for terminating the thermal pulse when the 

heating means 1 2 reaches a preselected maximum 
temperature; and temperature control means oper- 
ative^ connected to the heating means 12 and 
comprising temperature detecting means for de- 
io tecting the temperature of the heating means 12 
and feed back control means for controlling the 
heating means based on the detected temperature 
so that the temperature of the heating means is 
maintained in a preselected range. 

15 

6. A device as claimed in any one of claims 1 to 4, in 
which at least one sensor is positioned off-axis with 
respect to the longitudinal axis of the flow of fluid for 
detecting flow velocity along a line traversing the di- 

20 rection of flow of the fluid. 

7. A device as claimed in any of claims 1 to 6, in which 
two signal generating means 12 are provided, one 
on each side of the sensor(s) 14,16. 

25 

8. A device as claimed in any of claims 1 to 7, in which 
the detectable signal generated by the signal gen- 
erating means 12 comprises simultaneously at 
least one high frequency component and at least 

30 one lower frequency component. 

9. A device as claimed in any one of claims 4 to 8, in 
which the heating means 12 includes a source of 
power 10 for generating alternating power pulses to 

35 produce a detectable thermal signal. 

10. A device for measuring the flow of a fluid compris- 
ing: 

40 (a) detectable signal generation means for gen- 

erating a detectable signal and which inherently 
generates a crosstalk signal at half the frequen- 
cy of the detectable signal; 

45 (b) means for detecting the detectable signal 

downstream of the detectable signal genera- 
tion means; 

(c) signal process means for process the de- 
so tectable signal and obtaining a flow character- 
istic of the fluid therefrom; 

(d) crosstalk interference signal detection 
means for detecting the crosstalk interference 

55 signal generated by the detectable signal gen- 

eration means; and 

(e) means for eliminating any phase delay er- 
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rors between the detectable signals due to er- 
rors created by the components of the signal 
process means. 

1 1 . A device as claimed in any one of claims 1 to 1 0, in s 
which the or each signal generating means and/or 
the or each sensor comprises a thin film conductor 

32 comprising a solid electrically non-conducting 
core 34 and a thin metallic film 36 covering at least 
a portion of the core 34, the conductor being oper- 10 
atively connected to a source of electrical energy. 

12. A method of determining the time delay difference 
between two instances of the same signal unam- 
biguously, said signal containing at least two known is 
frequency components comprising: 

a) cross correlating the signal over a correlation 
interval substantially equal to the period of the 
lowest component; and 20 

b) locating the peak of the cross correlation 
function to thereby obtain the time delay. 
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